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We study inhomogeneous cosmological variations in the fine structure 'constant', a, in Friedmann 
universes. Inhomogeneous motions of the scalar field driving changes in a display spatial oscillations 
that decrease in amplitude with increasing time. The inhomogeneous evolution quickly approaches 
that found for exact Friedmann universes. We prove a theorem to show that oscillations of a in 
£N) ' time (or redshift) cannot occur in Friedmann universes in the BSBM theories considered here. 
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Elsewhere, |0,0|, we have discussed the behaviour of a class of cosmologies in an exact theory in which the fine 
structure "constant" varies in time. This theory of Sandvik, Barrow and Magueijo is an extension, to include the 
self-gravitation of the dielectric medium, of Bekenstein's prescription || for generalising Maxwell's equations to 
incorporate varying electron charge. We will refer to it as the BSBM theory. The fine structure "constant" a varies 
^v^j ' through the space-time dynamics of a scalar '"dielectric" field tp, (where a = exp[2?/>]) in these theories. However the 
^-H i overall behaviour is significantly affected by the form of the coupling. Even though the requirement that the energy in 
i\) be positive definite fixes the sign of the coupling constant uj we find that ip is driven by a term of the form £ em /uj, 
where C em is the electromagnetic Lagrangian. In general, C em can be positive or negative, a fact we parameterize in 
terms of ( = C e m/p, where p is the energy density. The sign of £ for the dark matter in the universe turns out to be 
\ of exceptional significance. 

■ In our earlier studies we have focussed on the case where ( < and the dark matter in the universe is dominated 
Oh' by magnetic field couplings. This was motivated by the discovery that C < matter leads to a slow (logarithmic) 
increase in the value of a with cosmic time during the matter era of the universe and constant behaviour during any 
period in which the expansion is curvature dominated, accelerates, or is dominated by radiation in universes with a 
matter-radiation balance like our own. Thus for £ < we find slow time-evolution of the fine structure 'constant' that 
is consistent with the observations of Webb et al of Aa/a = (-0.72 ± 0.18) x 10 -5 at z = 1 - 3.5. They are also 
consistent with the low-redshift, z ~ 0.1, upper limits on time variation of a provided by Oklo 0^],and high-redshift 
constraints imposed by the microwave background temperature fluctuations |],[l(|, and primordial nucleosynthesis 
p| . Other hints of varying constants in astronomical studies have recently been reported by Ivanchik et al |l^] . 

All of the studies described above have been performed in the context of an exact isotropic and homogeneous 
Friedmann universe. All variations in the fine structure "constant" therefore depend only on cosmic time. However, 
the rate of variation that is suggested by recent astronomical observations of quasar spectra, or allowed by geophysical 
data at recent times, is very small, Aa/a ~ 10 -5 , and spatial variations in the rate of time variation could easily be of 
similar order fl3|| . It is therefore important to determine if spatial variations in the rate of change of a are significant 
in the BSBM theory and whether they allow different modes of time variation to occur in addition to those studied 
in the purely homogeneous variations found in refs. [§,[[). 
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II. THE BSBM THEORY 



There are a number of possible theories allowing for the variation of the fine structure constant, a. In the simplest 
cases one takes c and h to be constants (see howev er p^| Jl6| ) and attributes variations in a to changes in the electron 
charge, e, or the permittivity of free space (see ref. KM for a discussion of the meaning of this choice). This is done by 
letting e take on the value of a real scalar field which varies in space and time (for more complicated cases, resorting 
to complex fields undergoing spontaneous symmetry breaking, see the case of fast tracks discussed in [|15[| ). Thus 
cq — > e = eoe(x 11 ), where e is a dimensionless scalar field and eo is a constant denoting the present value of e. This 
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operation implies that some well established assumptions, like charge conservation, must give way Nevertheless, 
the principles of local gauge invariance and causality are maintained, as is the scale invariance of the e field (under a 
suitable choice of dynamics). In addition there is no conflict with local Lorentz invariance or covariance. With this 
set up in mind, the dynamics of our theory is then constructed as follows. Since e is the electromagnetic coupling, 
the e field couples to the gauge field as eA^ in the Lagrangian and the gauge transformation which leaves the action 
invariant is eA^ — > eA^ + x,m> rather than the usual A^ — > A^ + x,m- The gauge- invariant electromagnetic field tensor 
is therefore 

= \ - M„),„) , (i) 

which reduces to the usual form when e is constant. The electromagnetic part of the action is still 

Sera = ~ [ d^X^F^F^, (2) 
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and the dynamics of the e field are controlled by the kinetic term 

(3) 

as in dilaton theories. Here, I is the characteristic length scale of the theory, introduced for dimensional reasons. This 
constant length scale gives the scale down to which the electric field around a point charge is accurately of Coulomb 
type. The corresponding energy scale, u> — Kc/l, has to lie between a few tens of MeV and Planck scale, ~ 10 19 GeV, 
to avoid conflict with experiment. 

Our generalisation of the scalar theory proposed by Bekenstein || described in ref. jij includes the gravitational 
effects of ip = loge. It gives the field equations: 

G, v = 8nG (T™ + T$ v + T^e" 2 *) , (4) 

where the various T M „ are the matter, ip and electromagnetic stress energy tensors. Recall the ip lagrangian is 
C^j, = — jd^ipd^ip an d the ip field obeys the equation of motion 

□V = -e-^em, (5) 
id 

where we have defined the coupling constant u> = (he) /I 2 . This constant is of order ~ 1 if, as in Q], the energy scale 
is similar to Planck scale. It is clear that C em vanishes for a sea of pure radiation since then C em — (E 2 — B 2 )/2 = 0. 
We therefore expect the variation in a to be driven by electrostatic and magnetostatic energy-components rather than 
electromagnetic radiation. In order to make quantitative predictions we need to know how much of the non-relativistic 
matter contributes to the RHS of Eqn. (||). This is parametrised by ( = £ em /p, where p is the energy density, and 
for baryonic matter C em = E 2 /2. For protons and neutrons Cp and Cn can be estimated from the electromagnetic 
corrections to the nucleon mass, 0.63 MeV and —0.13 MeV, respectively JT^JliJ, . This correction contains the E 2 /2 
contribution (always positive), but also terms of the form a M (where is the quarks' current) and so cannot be 
used directly. Hence we take a guiding value ( p w ( n ~ 10 -4 . Furthermore the cosmological value of £ (denoted 
Cm) has to be weighted by the fraction of matter that is non-baryonic, a point ignored in the literature (|J2^]. Hence, 
Cm depends strongly on the nature of the dark matter and can take both positive and negative values depending 
on which of Coulomb-energy or magnetostatic energy dominates the dark matter of the Universe. It could be that 
C cdm ~ —1 (superconducting cosmic strings, for which C em w — B 2 /2), or (cdm *C 1 (neutrinos). BBN predicts an 
approximate value for the baryon density of w 0.03 with a Hubble parameter of fig w 0.6 , implying Qcdm ~ 0.3. 
Thus depending on the nature of the dark matter Cm can be virtually anything between —1 and +1. The uncertainties 
in the underlying quark physics and especially the constituents of the dark matter make it difficult to impose more 
certain bounds on Cm- 

We should not confuse this theory with other similar variations. Bekenstein's theory |jj does not take into account 
the stress energy tensor of the dielectric field in Einstein's equations. Dilaton theories predict a global coupling 
between the scalar and all other matter fields. As a result they predict variations in other constants of nature, 
and also a different dynamics to all the matter coupled to clcctromagnetism. This model may be seen as a more 
conservative alternative to varying-speed-of-light scenarios p2],pT , p3| -p6| . An interesting application of our approach 



has also recently been made to braneworld cosmology by Youm |2 



Assuming a homogeneous and isotropic Friedmann 



metric with expansion scale factor a(t) and curvature parameter k in eqn. (|4j), we obtain the field equations (c = 1 
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( p m (1+ I Cm | exp [-2ip]) + p r exp [-2ip] + -ip' 
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where A is the cosmological constant. For the scalar field we have the propagation equation, 
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%j) + 3Hijj = exp [-2-ip]( m p m , (7) 

UJ 

where H = a /a is the Hubble expansion rate. Note that the sign of the evolution of ip is dependent on the sign of 
£ m . Since the observational data is consistent with a smaller value of a in the past, we will in this paper confine 
our study to negative values of £ m , in line with our recent discussion in Ref. . The conservation equations for the 
non-interacting radiation and matter densities are 

Pm + 3Hp m = (8) 

p r + AHp r = 2lpp r . (9) 

and so p m cx a~ 3 and p r e~ 2 ^ cx a~ 4 , respectively. If additional non- interacting perfect fluids satisfying equation of 
state p = (7 — l)p are added to the universe then they contribute density terms p cx a -37 to the RHS of eqn. (Q) as 
usual. 



III. INHOMOGENEOUS SOLUTIONS WITH VARYING a 



The Friedmann models with varying a have shown that when £ m < the homogeneous motion of the tp does not in 
general create significant metric perturbations at late times and we can safely assume that the expansion scale factor is 
that of the usual Friedmann universe for the appropriate fluid. The behaviour of ip then follows from a solution of the 
ip conservation equation in which the expansion scale factor is taken to be that of the Friedmann universe for matter 
with the same equation of state in general relativity {ip = £ = 0) . Our earlier analyses found that ip is approximately 
constant during the radiation era, and a increases as 2N hi(t) during the dust dominated era when spatial curvature 
is negligible, and tends to a constant in any subsequent era dominated by negative spatial curvature or a positive 
cosmological constant When ( < we can use the same test-motion approach to investigate inhomogeneous 
variations in ip and a as the universe expands. 

We assume that the expansion scale factor is that of the Friedmann model 

a = t n (10) 
and solve the wave equation in one of its appropriate forms: 

(11) 
(12) 
(13) 
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/3 m exp 
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0< 


) - aV 2 V> = 


ATexp[-2V>] 





where N is a constant, defined by 

N = -^ Pm a 3 > 0. 

UJ 

We can see in a general way that the effects of small inhomogeneities in the density of electromagnetically-coupled 
matter, will create a spatial variation in N(x) and this will create small spatial variations in a ~ 2N(x) \nt during 
the dust era. It is also possible for inhomogeneity in N to be created by variations in the value of £ for the form of 
matter that dominates on any particular scale. Our assumption of £ m < applies only to the dominant dark matter. 
On small scales the luminous matter might dominate and there will be a variation in the effective value (and even the 
sign) of C but we will not explore these possibilities further here. 

We seek a general solution of ( |l3| ) of the form 
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ip = iph + S(x,t) 



(14) 



where iph(t) is the solution to the space-independent problem (V^ = 0), so by definition iph(t) is an exact solution of 

(V^a 3 ) = Nexp[-2if> h ] 



d_ 
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We note immediately an important general property of this equation, that applies to all Friedmann universes with 
varying a : 

No- os dilation theorem: In the BSBM theory, a cannot display oscillatory behaviour in time in a Friedmann 
universe of any curvature. 

The proof is simple: When N is positive (negative) the right-hand side of equation (^|) is positive (negative), 
iph cannot have an expansion maximum (minimum) since iph — and iph < (> 0) there. Therefore iph cannot 
oscillate in time and so neither can a — exp[2ip]. 

We see that in the case of interest, when N > 0, tp can have a minimum but thereafter it must always increase 
irrespective of the behaviour of the expansion scale factor. However, if the equation is linearised in iph this is no longer 
true if attention is not confined to the small ip regime where exp[— 2iph] ~ 1 — 2tph > and spurious oscillations of 
ip (and a) in time can appear to arise at late times if ip grows. It is of particular interest that this proof that iph 
cannot have a maximum applies to recollapsing universes (k = +1) as well as to ever-expanding universes {k < 0). It 
also means that oscillations of a with redshift should not be observed in Friedmann universe. This might prove an 
interesting prediction for future observations to test. 

Substituting Q into (jT|) we get 



dt 



((5a 3 ) - aV 2 <5 = N exp[-2t/> h ]{exp[-2<J] - 1} 
So for small 5 

— (<5a 3 ) - aV 2 5 = -2N5exp[-2ip h ] + 0(5 2 ) 

Now look for separable solutions 



5 = T(t)D(x) 



and we have 



T T 2N V D 

a 2 + 3aa- + exp[-2^] - -fi 2 - — — (15) 

11a D 

where /i 2 is a separation constant with a sign chosen to ensure non-growing, oscillatory, inhomogeneity in D{x) at 
spatial infinity . In this equation we can always neglect 2Na~ l exp[— 2iph] with respect to /i 2 as t — * oo because iph 
never falls with time (in the dust era iph grows as \ In [2 AT ln(£)] as t — > oo, for example, This is an important 
feature of the variation of ip, and a, in BSBM varying-a theories when £ < 0. It ensures that the kinetic term and 
the Cmexp[— 2ip] terms can be neglected in the Friedmann equation asymptotically and the expansion scale factor 
can self-consistently be assumed to be of the same form as when a does not vary (this is not true if £ > 0). Note 
that in the Friedmann case (5 = 0) we can evaluate the corrections to the test-motion approximation by calculating 
the leading order corrections to the Friedmann equation if we use the solution for ip found from the solution of the 
wave equation. These are largest for the dust universe but even there we find that the next order correction to the 
first-order assumption that a(t) = t 2 / 3 is a(t) = i 2/3 (lni) with |C| /3 - 0.3 - 0.03 and so are small. 
Hence, in this approximation we have 

•■ 3a • u 2 T , s 

T+—T+^- = Q 16 
a a z 

and 

V 2 L> = -fi 2 D 

so we have the standard separable spherical oscillator solution in spherical polar coordinates: 
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D(r,0,<p) = Y. c ^ Y c( 9 ^)r- 1/2 Z e+ i(Lir) 



where Z is a cylindrical function and Y the spherical harmonic function. If we specialise to spatially-flat cosmologies 
with perfect fluid equations of state for pressure p and density p of the form 

P = (7 - !)a 

then the expansion scale factor will have power law form 

a = t n (17) 

with n = 2/37. In these cases we have 

tf + 3nT + /^Ti 1 " 2 " = (18) 

We are interested in inhomogeneous solutions which introduce new behaviour as a result of including inhomogeneity. 
Therefore we impose a boundary condition that T = for /1 = since when D = the time-dependent solution is 
already included in iphit). Thus, for n < 1 : 



T(t) = tV-*"V*C x Z v ( -^t 1 -' 
\ 1 — n 

_ |l-3n| 



(19) 
(20) 



2(1 -n) 

where Z(..) is a cylindrical function, while for the curvature-dominated expansion with n = 1 : 

Toe* 9 (21) 
9= -1± a/1 -M 2 (22) 

and we choose the + solution to satisfy the boundary condition. The late-time behaviour is easily determined as 
t — > 00: 

T(t) (x r n x oscillations; n ^ 1. (23) 



TO) oc t-i+\/w^ ; n = 1 (24) 



and decays, T oc a -1 , as £ — * 00. However, as we have already pointed out the oscillatory behaviour is an artefact of 
the linearisation process and the Bessel-like oscillation are not reached by the solution for ip. 

In the radiation era we can find an solution of eqn. ( p"5| ) for T(t) without neglecting the term 2Na~ 1 expj— 2iph] 
since the radiation universe has the simple exact solution: 

4> h = ± log(8JV) + \ log(t) (25) 

Substituting ( p5| ) in eqn. ( p^| ) we find 

T W = ^{^m(2^ 1/2 )+5J_ m (2 M t 1 / 2 )} 

where 

iV3 

m = 

2 

and we see explicitly that there is agreement with the asymptote ( |23| ) of the approximated equation when n = 1/2. 
The boundary condition for transition to the homogeneous problems requires that we put B = and again the late- 
time oscillations are recognised as arising purely from the linearisation process. Similar exact solutions can be found 
for all universes with 1/3 < n < 2/3. 

The cosmological constant case of 7 = is distinct, with 
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a — cxp[Hot] 



which gives 

= f + 3HqT + /i 2 T exp[-2H t] ~ f + 3H f 

as t — > oo, so 

T -> ^ - rrr exp[-3JIo(i + io)] -» 

MO 

This behaviour is in accord with the expectations of a cosmic no hair theorem. It means that if a period of inflation 
occurs in the very early universe then large scale inhomogcncity will appear increasingly negligible with time within 
the event horizon of a geodesically moving observer. In the late stages of a universe like our own, which displays 
evidence of being accelerated by the presence of a positive cosmological constant, (27), it ensures that time variations 
in a will not grow. This is to be expected since the inhomogeneities in density are also prevented from growing by 
the effects of the cosmological constant. 



IV. THE CASE OF C > 

When the dark matter is dominated by electric field energy, we have £ > 0, and the behaviour of eq. (fjj) is 
very different to that obtained when £ < 0. Most crucially, the test-motion approximation used above to analyse 
the behaviour or (0) does not apply, even for the purely time-dependent ip evolution in a Friedmann universe. The 
solutions obtained for ip by assuming the scale factor evolution a = t n of general relativity (with constant a) lead 
to solutions for ip (and a) which do not increase with time. For example, we have a a T 1 in the curvature era 
and a oc ln(ioA) in the dust era. These contribute kinetic {ip 2 ) and magnetic contributions (Cexp[— 2ip}) terms 
which dominate the underlying Friedmann equation, (^), at large times and the expansion of the universe is not well 
approximated by that obtained in general relativistic cosmologies with the same equation of state and constant a 
except over finite non-asymptotic intervals of time. This leads to problems accommodating observational constraints, 
notably the results of studies of the structure of the microwave background at last scattering 

§H and big bang 



nucleosynthesis ]ll|] in the radiation era because the value of a then is significantly different from today, unlike in the 
cases of ( < 0. Cosmologies with £ > have been discussed in ref. jl9) in a theory that is similar in structure to the 
BSBM theory discussed here. We will discuss the ( > version of the theory in more detail elsewhere. It is less well 
behaved and does not seem to provide the smooth and simple perturbation of the standard cosmology with constant 
a as seen in the negative C case. 



V. DISCUSSION 



We have shown that the time-dependent solutions to the Friedmann model are stable against the effects of inhomo- 
geneous motions of the ip field. In the case of inhomogeneous variation the cosmological solutions in universes with 
scale factor ait) = t n , to leading order take the form: 

lP(x,t) = Tjj h (t) + C^ 1 - 3 "^ 2 

when n ^ 1, and 

iP(x,t) = ip h (t) + At~ 1+ V^ ^^(^K'/^^ifH 

1=0 

when n = l, while for the case of a = exp[iJoi] : 

ip(x,t) = ip h {t) + O(exp[-3H t ]). 

Thus in all cases we have 



Ju 



1 - n 



E 

n n 



Apr) 



tp(x,t) — > ip h (t) 



G 



as t — > oo and at late times spatial variations in the fine structure constant decay as 

oo 

a = cxp[2^]{l + 2r n Y / c^eYe{0,tp)r- 1/2 Z e+ i( f ir) + ..} 
for n =/= 1. Hence, denoting ah = exp[2i/jft], the spatial variation in a decays in time in the n^l universes as 

a a/* 2 

Analogous expressions can be written down mutatis mutandis for 8a /a in the n = 0, 1 cases. 

It is important to compare the evolution of the fine structure constant a(t) in the BSBM theory in the homogeneous 
case with that for the situation admitting inhomogeneous motions of the fine structure 'constant', a(t,x), here. 
To leading order, the overall pattern of time evolution studied in refs. |j]J^] is unaffected by the presence of small 
inhomogeneities. However, small spatial variations of an oscillatory character are expected to exist in the value of 
the fine structure constant over astronomical scales, reflecting the non-linear self-interaction of the a{ip) field which 
carries the variations in a. The spatial variation amplitudes, 8a/ a, are found to decay with time as the universe 
expands and will not be as significant as the overall variation in time of the mean value of a(t) oc ln(t) during the 
dust-dominated phase of a spatially-flat universe. Inhomogeneous test motions of the ip field will have been decaying 
in amplitude throughout the period when the universe was dominated by dust if £ < 0. Therefore we would not expect 
any significant inhomogeneities to survive at the astronomically interesting epoch z ~ 1 — 4 where the value of the fine 
structure constant can be probed spectroscopically with high precision. However, our discussion has not considered 
three situations where more significant spatial variations might arise. The first is the situation within gravitationally 
bound matter inhomogeneities of large scale which separate out from the expansion of the Universe and collapse to 
form superclusters and clusters of galaxies. These behave in a manner similar to that expected of separate closed 
universes until deviations from spherical symmetry become significant. Our analysis is not applicable here because the 
dynamics of the bound inhomogeneities will differ significantly after they separate off from the background expansion. 
The second situation of interest is that in which perturbations of the Friedmann metric are included in the problem 
and allowed to couple to spatial variations in ip, or a. This coupling will lead to small temperature fluctuations in the 
microwave background radiation. Finally, the variation in the value and sign of £ with scale for the dominant form of 
matter could introduce a distinctive inhomogeneity. These problems will be discussed elsewhere. 
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